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Outflows from Quasars and Ultraluminous X-ray Sources 

A.R. King a 
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Mass outflows from Eddington-limited accreting compact objects appear to be a very widespread phenomenon. 
They may provide the soft excess observed in quasars and ULXs, and imply that such objects have a major effect 
on their surroundings. In particular they allow a simple parameter-free argument for the Mbh — o relation for 
galaxies, and offer a straightforward interpretation of the emission nebulae seen around ULXs. 



1. INTRODUCTION 

Recent XMM-Newton observations of bright 
quasars (Pounds et al., 2003a, b; Reeves et 
al., 2003) give strong evidence for intense out- 
flows from the nucleus, with mass rates M out ~ 
lM© yr^ 1 and velocity v ~ 0.1c, in the form of 
blueshifted X-ray absorption lines. These out- 
flows closely resemble those recently inferred in a 
set of ultraluminous X-ray sources (ULXs) with 
extremely soft spectral components (Mukai et al, 
2003; Fabbiano et al., 2003). Simple theory shows 
that the outflows are probably optically thick to 
electron scattering, with a photosphere of ~ 100 
Schwarzschild radii, and driven by continuum ra- 
diation pressure. In all cases the outflow velocity 
is close to the escape velocity from the scatter- 
ing photosphere. As a result the outflow momen- 
tum flux is comparable to that in the Eddington- 
limited radiation field, i.e. 



M out v 



L 



Edd 



(1) 



where M out is the mass outflow rate and LEdd 
the Eddington luminosity, while the mechanical 
energy flux is 



^Moutv 2 



^Edd 

2M out c 2 ' 



(2) 



These outflows appear to be a widespread phe- 
nomenon, not only in currently observed sys- 
tems such as quasars and ULXs, but also in the 
growth of supermassive black holes in the centres 
of galaxies in the past. 



2. OUTFLOWS FROM EDDINGTON 
LIMITED ACCRETORS 

I outline here a simple theory of outflows from 
black holes accreting at rates comparable with the 
Eddington value 



M, 



Edd 



4ttGM 

7JKC 



(3) 



Here rjc 2 is the accretion yield from unit mass, 
and k is the electron scattering opacity. We as- 
sume that the outflow is radial, in a double cone 
occupying solid angle Airb, and has constant speed 
v for sufficiently large radial distance r. I will jus- 
tify the second assumption later in this Section. 
Mass conservation implies an outflow density 



Mo 



Anvbr 2 



(4) 



The electron scattering optical depth through 
the outflow, viewed from infinity down to radius 
R, is 



nM out 

npar = — . 

R r invbR 

From we get 

- 1 c M out 
2r)b R v A/ Edd ' 



(5) 



(6) 



Defining the photosphcric radius R p h as the point 
r = 1 gives 
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P h 



1 c M n 



5 c M 



R s 2rjb v M E dd & M Ed d 



(7) 



2 



where we have taken r\ ~ 0.1 at the last step. If 
b < l,v/c < 1 we see that i? p h > R s for any 
outflow rate M out of order M Edd , that is, such 
outflows are Compton thick. If instead b « 1, 
photons typically escape from the side of the out- 
flow rather than making their way radially out- 
wards through all of it. Almost all of the photons 
escape in this way within radial distance r = R± 
where the optical depth across the flow 



T_|_ 



np{r)b 



1/2, 



is of order unity. Thus 

R±_ _ 1 c Mout ^ 
R s ~ 2776 1 / 2 v M Edd " 
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(8) 
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and we again conclude that the outflow is Comp- 
ton thick for M out ~ M Ed d- 

This evidently implies that much of the emis- 
sion from such objects will be thermalized and 
observed as a softer spectral component (see eq 
EOlbelow) . The observed harder X-rays must pre- 
sumably either be produced near the skin of the 
outflow (i.e. at moderate r), or result from shocks 
within the outflow. In both cases their total lumi- 
nosity must be lower than that of the thermalized 
soft component. 

We now investigate how the outflow is driven. 
Since the wind is Compton thick most of the pho- 
tons have scattered and thus on average given up 
their original momentum to the outflow. Out- 
side the radius i? p h or R± the photons decouple 
from the matter and there is no more accelera- 
tion. This justifies our assumption that v is con- 
stant for large r, and it is self-consistent to use 
the assumption to integrate inwards to i? p h or 
R ± . 

To ensure that the matter reaches the escape 
speed we assume that the radii R p h,R± to lie 
close to the escape radius i? csc , i-e. 



-Rph,-R_L — Resc = ~^Rs- 



(10) 



If instead the photosphere is outside the escape 
radius, i.e. there is substantial photon trapping, 
one would insert an optical depth factor r on the 
rhs. The quasar observations referred to above, 
and the normalization of the A/bh — c relation 



derived below both suggest that r ~ 1, and I 
adopt this value in the rest of this article. 
From JTHJ and EI) we find 
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M out V27?/M Edd 

where / = 6, b 1 / 2 in the two cases b < 1, b << 1 
We can write these formulae more compactly as 

v _ 2/L Edd 

c ~ MoutC 2 
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(12) 
(13) 



We note that / ~ 1 except for very narrowly col- 
limated outflows (b < 10~ 2 ). 

An immediate consequence of <|12|) is 



M out v = 2/ 



^Edd 



(14) 



i.e. the momentum flux in the wind is always of 
the same order as that in the Eddington-limited 
radiation field, as expected for an Compton thick 
wind driven by radiation pressure. The energy 
flux (mechanical luminosity) of the wind is lower 
than that of the radiation field by a factor of order 
v/c: 



v. T 2(/L Ed d) s 

-/^Edd - — ~ — 

c M out c 2 



(15) 



3. THE BLACKBODY COMPONENT 

Since the outflow is Compton thick for M ou t ~ 
M Edd , much of the accretion luminosity generated 
deep in the potential well near R s must emerge 
as blackbody-like emission from it. If b ~ 1 the 
quasi-spherical radiating area is 



.4 



phot 



= AirbR. 



ph 



(16) 



If instead b << 1 the accretion luminosity 
emerges from the curved sides of the outflow 
cones, with radiating area 



A ± = 2irR ph 



> 6 V2 



We can combine these two cases as 
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with g(b) = 1/6, l/2fo 1//2 in the two cases. Again 
<7 ~ 1 unless b < 10~ 2 , so the areas are similar 
in the two cases. However we note that the ra- 
diation patterns differ. In particular if b is small 
radiation is emitted over a wider solid angle than 
the matter. Numerically we have 

A ph<1 _ = 2 x 10 29 gM*M- 2 cm 2 , (19) 

where M x = M out /(1M yr^ 1 ), M 8 = M/1O 8 M . 
The effective blackbody temperature is 

T cff = 1 x 10V 1/4 ^T lM $ /4 k - (20) 

Clearly such a component is a promising candi- 
date for the soft excess observed in many AGN 
and ULXs. 

4. THE M BH - (J RELATION 

It is now widely accepted that the centre of 
every galaxy contains a supermassive black hole. 
The close observational correlation (Ferrarese & 
Merritt, 2000; Gebhardt et al., 2000; Tremaine et 
al., 2002) between the mass M of this hole and the 
velocity dispersion a of the host bulge strongly 
suggests a connection between the formation of 
the black hole and of the galaxy itself. 

The outflows considered above must have an 
important effect on this problem. We know that 
most of the mass of the nuclear black holes is as- 
sembled by luminous accretion (Soltan 1982; Yu 
& Tremaine, 2002). It seems likely that the rate 
at which mass tries to flow in towards the central 
black hole in a galaxy is set by conditions far from 
the hole, for example by interactions or mergers 
with other galaxies. It is quite possible there- 
fore that super-Eddington conditions prevail for 
most of the time that the central black hole mass 
is being built up. Evidently the resulting Edding- 
ton thrust l|14|) can have an important effect on 
the host galaxy. Unlike luminous energy, a large 
fraction of a mechanical energy flux like (|15[) is 
likely to be absorbed within the galaxy. To reach 
its present mass the black hole in PG1211+143 
could have accreted at a rate comparable to its 
current one for ~ 5 x 10 7 yr. During that time, 
an outflow like the observed one could have de- 
posited almost 10 60 erg in the host galaxy. This 



exceeds the binding energy ~ 10 59 erg of a bulge 
with 10 11 M and a ~ 300 km s" 1 . 

Ideas presented by Silk and Rees (1998, hence- 
forth SR98) and also considered by Haehnelt et 
al. (1998), Blandford (1999) and Fabian (1999) 
are relevant here. These authors envisage a situ- 
ation in which the initial black holes formed with 
masses ~ 1O 6 M before most of the stars. Ac- 
cretion on to these black holes is assumed to pro- 
duce outflow, which interacts with the surround- 
ing gas. Without a detailed treatment of the out- 
flow from a supercritically accreting black hole, 
SR98 used dimensional arguments to suggest a 
relation between M and a. However this still has 
a free parameter. Given the simple relation i|14fl 
one can now remove this freedom. The situation 
turns out to resemble a momentum-driven stellar 
wind bubble. Modelling this gives an Mbh — f re- 
lation devoid of free parameters, and remarkably 
close to the observed relation. 

4.1. Black hole wind bubbles 

SR98 modelled a protogalaxy as an isothermal 
sphere of dark matter. If the gas fraction is f g = 
^baryon /^matter ^ 0.16 (Spergel et al., 2003) its 
density is 



where a is assumed constant. The gas mass inside 
radius R is 

M(R) = ^ [ R pr 2 dr = (22) 
Jo " 

I now assume that mass flows towards the central 
black hole at some supercritical rate M acc , and 
thus exerts a momentum flux l|14|) on the sur- 
rounding gas, sweeping it up into a shell. As is 
well known from the theory of stellar wind bub- 
bles (e.g. Lamers & Casinelli 1999) this shell is 
bounded by an inner shock where the wind ve- 
locity is thermalized, and an outer shock where 
the surrounding gas is heated and compressed by 
the wind. These two regions are separated by 
a contact discontinuity. The shell velocity de- 
pends on whether the shocked wind gas is able to 
cool ('momentum-driven' flow) or not ('energy- 
driven' flow). In the absence of a detailed treat- 
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ment of a quasar wind, SR98 appear to have as- 
sumed the second case. In fact for the supercriti- 
cal outflows envisaged here it is easy to show that 
the bubble is efficiently Compton cooled, and thus 
momentum-driven instead (King, 2003) 

4.2. The M B h cr relation 

The speed v m of the momentum-driven shell 
now follows from the standard wind bubble argu- 
ment. At sufficiently large radii R the swept- 
up shell mass M(R) is much larger than the 
wind mass, and the shell expands under the im- 
pinging wind ram pressure pv 2 (this character- 
izes momentum-driven flows; in an energy-driven 
flow the thermal pressure of the shocked wind gas 
is dominant, while in a supernova blast wave the 
momentum injection is instantaneous rather than 
continuous). The shell's equation of motion is 
thus 



d 
dt 



M(R)R 



4irR 2 pv 2 = M a 



L 



Edd 



(23) 



where we have used first the mass conservation 
equation for the quasar wind, and then (|14() to 
simplify the rhs. Integrating this equation for R 
with the final form of the rhs gives 



M (R)R = 



(24) 



where I have neglected the integration constant 
as M(R) is dominated by swept-up mass at large 
t. Using (1221) for M (R) and integrating once more 
gives 



R A 



G^Edd 2 



(25) 



where again we may neglect the integration con- 
stant for large t. We see that in the snow- 
plow phase the shell moves with constant velocity 
v m — R/t, with 

2 GiEdd 



2/»oV 



(26) 



Note that this velocity is larger for higher LEdd, 
i.e. higher black hole mass. This solution holds if 
the shell is inside the cooling radius R c ; outside 
this radius the shell speed eventually increases to 
the energy-driven value v e , which also grows with 
M BH . 
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1.5 x 10V 



8„4 



200 
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(27) 



I now consider the growth of the black hole 
mass by accretion. Initially the mass is small, 
inflow is definitely supercritical, and even the 
energy-driven shell velocity would be smaller 
than the escape velocity a. No mass is driven 
away, and accretion at a rate MEdd can occur effi- 
ciently. However as the black hole grows, we even- 
tually reach a situation in which v e > a > v m . 
Further growth is now only possible until the shell 
reaches R c , and then only until the point where 
Vm = c. Thus given an adequate mass supply, 
e.g. through mergers, the final black hole mass 
is given by setting v m = a in (|26|l . Thus wc find 
the relation 

fg K 4 

This is remarkably close to the observed relation 
(Tremaine et al., 2000). The relation also has no 
free parameter, unlike earlier derivations (SR98; 
Haehnelt et al., 1998; Blandford 1999; Fabian 
1999) If v had been larger by an optical depth fac- 
tor t > 1 (i.e. i? p h, R± were outside the escape 
radius -Rose) a factor 1/r would have appeared 
on the rhs. Taking this factor ~ 1 as implic- 
itly done in deriving eqn llOfl above agrees with 
the observed outflow velocities ~ {GM / R^) 1 / 2 
in supercritically accreting quasars. 

This derivation of the Mbh — & relation implies 
that the central black holes in galaxies gain most 
of their mass in phases of super-Eddington in- 
flow. As relatively few AGN are observed in such 
phases, these must either be obscured (cf Fabian, 
1999) or at high redshift. It appears then that 
those quasars which are apparently now accreting 
at such rates (Pounds et al., 2003a, b; Reeves et 
al., 2003) are laggards in gaining mass. This idea 
agrees with the general picture that these objects 
- all narrow-line quasars - are super-Eddington 
because they have low black-hole masses, rather 
than unusually high mass inflow rates. 

5. ULX NEBULAE 

We have seen that outflows from Eddington- 
limited AGN can have a profound effect on their 
surroundings, in the form of the Mbh - f relation. 
In the stellar-mass context, ULXs are the most 
supercritically accreting objects, and one would 
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naturally expect a similar effect here too. Indeed, 
many ULXs are observed to lie in the centres 
of unusually large and bright emission nebulae. 
These have in the past been variously interpreted 
as supernova/hypernova remnants (cf Roberts et 
al, 2003) or photoionized nebulae (Pakull & Miri- 
oni, 2003). However it is clear from the work 
of the last Section that a more likely explana- 
tion of the unusual energetics of these objects is 
as bubbles blown in the ISM by an Eddington 
outflow from the ULX. Pakull & Mirioni (2003) 
indeed note that the nebulae are probably shock 
excited. This interpretation means that the neb- 
ulae are powered by the accumulated kinetic en- 
ergy of the outflow that was dissipated over the 
gas cooling time in the nebula (or the ULX life- 
time if this is shorter). In particular there is no 
requirement that the radiation currently emitted 
from the ULX should power the emission nebula. 
The nebulae do not provide arguments against 
anisotropic emission from this source: even if the 
wind is anisotropic, it is probable that it would 
have precessed many times over the ULX's life- 
time, and thus left a roughly spherical nebula, as 
observed. One cannot therefore use ULX nebu- 
lae to argue against the interpretation of ULXs as 
anisotropically emitting stellar-mass X-ray bina- 
ries (cf King et al., 2001). 



6. CONCLUSION 

Mass outflows from Eddington-limited ac- 
creting compact objects appear to be a very 
widespread phenomenon. They have major im- 
plications for quasars and galaxy formation, and 
for ULXs. Some of these are barely explored as 
yet, and the field promises to be fruitful. 
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